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ABSTRACT 



In a previous paper (Serna & Alimi 1996b), we have pointed out the existence of some 
particular scalar-tensor gravity theories able to relax the nucleosynthesis constraint on the cosmic 
baryonic density. In this paper, we present an exhaustive study of primordial nucleosynthesis 
in the framework of such theories taking into account the currently adopted observational 
constraints. We show that a wide class of them allows for a baryonic density very close to that 
needed for the universe closure. This class of theories converges soon enough towards General 
Relativity and, hence, is compatible with all solar-system and binary pulsar gravitational tests. 
In other words, we show that primordial nucleosynthesis does not always impose a very stringent 
bound on the baryon contribution to the density parameter. 

Subject headings: cosmology: early universe, cosmology: dark matter, cosmology: theory 
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1. INTRODUCTION 



One of the most convincing pieces supporting the standard hot big bang model for the early Universe is 
the excellent agreement between the predicted and observed primordial light element abundances. However, 
this agreement requires that the present baryon density pb must be smaller (Walker et al. 1991, Copi, 
Schramm & Turner 1995a-1995b, Olive 1996), than about 10% of that needed for the Universe closure 
(fig = Pb/pcrit < 0.1). A high-density universe (O 3> 0.1), seems then to be impossible unless that most 
dark matter in the Universe were non-baryonic. 

The problem in admitting large pb values is that, in the standard primordial nucleosynthesis model, 
they lead to the overproduction of 4 He and 7 Li and the underproduction of deuterium. In order to 
understand this difficulty, we recall that primordial nucleosynthesis starts, in the early Universe, soon after 
that the cosmic temperature becomes smaller than that needed to maintain the proton-to-neutron ratio 
in its equilibrium value (freezing-out temperature). The nuclear reactions which then take place lead first 
to the D and 3 He formation. These elements are then burnt to produce 4 He. At the end of this process, 
essentially all neutrons are incorporated into 4 He, while just a small fraction of them remains in D and 3 He. 
Consequently, the He abundance mainly depends on the fraction of neutrons existing at the beginning of 
the nucleosynthesis process or, equivalently, on the freezing-out temperature for the n/p ratio. On the other 
hand, since the burning of D and 3 He grows with the baryon density, the final yields of these two elements 
decrease with pb, while those of 4 He and 7 Li increase^. 

Several attempts to relax the nucleosynthesis bound on Qb have been previously considered (Malaney 
& Mathews 1993). This is the case, for example, of inhomogeneous big bang models, and the attempts of 
modifying the nuclear reaction rates by the introduction of new decaying particles (e.g. Schramm 1991, 
Gyuk & Turner 1994). However, all these attempts have been unsuccessful or seem to be mere modifications 
to try to solve this problem. The nucleosynthesis bound on Qb is today considered as a unavoidable 
constraint which must be imposed even to study the formation of large-scale structures in the Universe. 

Another possibility also analyzed in the literature consists of modifying the gravity description and, 
hence, the Universe expansion rate without altering the nuclear physics, the Universe composition or the 
conservation laws. This last possibility can be considered as more natural from a theoretical point of view 

1 Although the lithium-7 history is much more complicated, when the baryon density is relatively large, 
its final abundance also increases with pb due to the 7 Li production through 7 Be 
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because it would allow us to construct a global framework which is coherent with that often assumed to 
describe the epochs prior to primordial nucleosynthesis. As a matter of fact, to study the earliest times 
in the Universe evolution, the existence of an extra scalar field, in addition to the Einstein metric tensor, 
is usually assumed. Such gravity theories are termed scalar-tensor theories (Bergmann 1968, Wagoner 
1970, Nordtvedt 1970), and arc the simplest generalization of the Einstein theory. They provide a natural 
(non-fine-tuned) way to restore the original ideas of inflation while avoiding the cosmological difficulties 
coming from the vacuum-dominated exponential expansion obtained in General Relativity (GR) (La & 
Stcinhardt 1989, Weinberg 1989, Steinhardt & Accettta 1990, Barrow & Maeda 1990, Liddle & Wands 1992, 
Deruelle, Gundlach & Langlois 1992, Garcia-Bellido & Wands 1992, Barrow 1995). Scalar-tensor theories 
also arise in the current theoretical attempts at deepening the connection between gravitation and the other 
interactions. For example, in modern revivals of the Kaluza-Klein theory and in supersymmetric theories 
with extra dimensions, one or more scalar fields arise in the compactification of these extra dimensions 
(Jordan 1949, Kolb, Perry & Walker 1986, Cho 1992, Wesson & Ponce dc Leon 1995). Furthermore, 
scalar-tensor theories may also appear as a low-energy limit of superstring theories (Green, Schwarz & 
Witten 1988). 

The Big Bang nucleosynthesis (BBN) in presence of a scalar field has been previously analyzed in some 
cases (Bekenstein & Meiseils 1980, Arai, Hashimoto & Fukui 1987, Serna, Dominguez-Tenreiro & Yepes 
1992, Serna, Dominguez-Tenreiro 1993, Serna & Alimi 1996b). Among all these analyses, only some models 
pointed out by Serna & Alimi (1996b) were able to relax the nucleosynthesis bound on Qb- The key of such 
a result was that the expansion rate resulting in such theories avoided simultaneously the overproduction of 
4 He and 7 Li and the underproduction of D. In this paper we construct a wide class of scalar-tensor gravity 
theories also allowing a large range for baryonic density. 



2. SCALAR-TENSOR NUCLEOSYNTHESIS 



Scalar-tensor theories are characterized by an arbitrary coupling function w($), which determines 
the relative importance of the additional scalar field $. GR is the asymptotic case in which the coupling 
function is infinite and $ = constant = 1. Since GR reproduces very accurately the solar-system and 
binary pulsar dynamics (Will 1993), the behavior of any viable scalar-tensor theory must be restricted to 
be at present extremely close to that implied by GR. The physical conditions in the early Universe and 
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the subsequent cosmological evolution in the framework of these theories can be nevertheless very different 
from the usual ones. 

In order to analyze the possible existence of scalar-tensor theories, able to reproduce the right primordial 
abundances for an interval much wider than in GR, we must first calculate the resulting cosmological 
models. To that end, we consider an isotropic and homogeneous universe and a specific form for the 
coupling function. The line element has then a Robertson- Walker form and the energy-momentum tensor 
corresponds to that of a perfect fluid. In that concerning the coupling function, a convenient form is given 

by 

\'S + 2uj\= (S/A 2 )^- 1 +k) (1) 

where A and k are arbitrary constants and x =| <& — 1 1. As a matter of fact, as shown by Serna & Alimi 
(1996a), such a form gives an exact representation for most of the particular scalar- tensor theories proposed 
in the literature and, in addition, it contains all the possible early behaviors of any theory where w($) is a 
monotonic, but arbitrary, function of $. 

When such an w($) function is introduced into the field equations, one essentially obtains four different 
classes of scalar-tensor theories. The two first classes correspond to singular models with a monotonic time 
evolution of the speed-up factor (£ = H/H FRW , where H is the Hubble parameter, while H FRW is that 
predicted by GR at the same temperature). In the first class, the Universe expansion is always faster than 
in GR (£ > 1), while the second class is characterized by an expansion rate slower than in GR (£ < 1). 
Nonsingular models or models with a critical temperature where 3 + 2uj = constitute the third class 
of theories. Finally, the fourth class corresponds to models with a non-monotonic £(T) function. These 
last theories have an initial phase where the Universe expansion is slower than in GR but, afterwards, it 
becomes faster than in the standard cosmology. Obviously, the Universe evolution in the framework of any 
viable scalar-tensor theory must satisfy £ — > 1 at present. 

In the first three classes of scalar-tensor theories, the right primordial yields are only obtained (Serna 
& Alimi 1996b) in the limit where such theories are almost indistinguishable from GR, from the beginning 
of primordial nucleosynthesis up to the present. Consequently, the allowed Ob interval is essentially the 
same as that implied by GR. 

We will then focus on the primordial production of light elements in the framework of theories with a 
non-monotonic evolution of the speed-up factor (class-4 models). This class of models is inevitably obtained 
when X/Vk > 1, and both 3 + 2lo and ($ — 1) have positive values. We show below that new constraints on 
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the baryonic density parameter can then be obtained. 

The computation of the primordial production of light elements^] has been performed by using the 
updated reaction rates of Caughlan and Fowler (1988) and Smith et al (1993). We have considered 
N v = 3, To = 2.73K, r„ = 889s, Ho = 50km s^Mpc -1 . Other runs have been also performed, with 
r n = 889.8 ± 3.6s and Hq = 100km s _1 Mpc _1 , in order to test the robustness of our conclusions. For each 
theory defined by the parameters A and k, we have searched the largest f2p values for which there exists a 
present value of the coupling function ujq (i) compatible with all solar system experiments and (ii) leading 
to the observed primordial abundances (Copi, Schramm & Turner 1995b): 

0.221 < Y P < 0.243 
D/H > 1.5 • 10~ 5 
(D+ 3 He)/H < 1.1 • 10~ 4 
0.7 ■ 10~ 10 < 7 Li/H < 3.5 ■ 10- 10 

Some of these largest upper bounds on are shown as isocontour lines on figure 1. The luq values 
satisfying the two previous criteria are different for different theories and, in particular, they are smaller 
when the upper bound on f2p are larger (see table) . The baryon contribution to the total density parameter 
cannot be therefore arbitrarily high because, for too large f2p values, the resulting cosmological models will 
not be then compatible with solar system experiments, which require loq > 500. The allowed interval found 
for fip is 

0.01 <Q B < 0.78 

When a wider observational range for the 7 Li/H abundance is considered ( 7 Li/H < 6 • 10~ 10 , which 
corresponds to a depletion by a factor of 4), this interval becomes 

0.01 <n B < 0.98 

In order to illustrate more precisely our procedure, we plot on figures 2, the light elements primordial 
abundances for GR and for two particular theories appearing on figure 1 . These two theories are defined 

2 Our code solves the scalar-tensor cosmological equations by using a 6 th order Runge-Kutta integration, 
while the primordial nucleosynthesis rate equations are integrated by the Beaudet & Yahil (1977) scheme. 
The outputs of this code were extensively tested and compared to those obtained by other authors, as well 
as to the analytical solutions known for some particular cases. It was then used in previous works as, e.g., 
Serna et al. (1992), Serna & Dominguez-Tenreiro (1992, 1993), Serna & Alimi (1996a, 1996b) 
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by A 2 = 0.4, logio(A 2 /fc) = 7.8 and oj q = 1.3 • 10 9 (square symbol), A 2 = 0.2, log 10 (\ 2 /k) = 8. and 
ujo = 5.75 • 10 3 (circle symbol). On these figures we see clearly that a range of large Qb values exists, where 
predicted primordial abundances are all compatible with observations. Any smaller upper bound on Qb 
than those seen on these figures, can be however obtained by considering larger oj values. We note also 
that, in fact, the most constraining light elements are 4 He and 7 Li. Any modification on the observational 
limits on these two elements would imply a different allowed range for fls- In particular, we see from figure 
2d, that an almost flat Universe ~ 1) is permitted when the upper observational limit for 7 Li/H is 
6-Hr 10 . 



3. Discussion and Conclusions 

A homogeneous and isotropic universe with a baryonic density very close to that needed for closure 
by baryons is then possible in the framework of class-4 scalar-tensor theories even when the Universe is 
composed solely by known particles. 

Since ojq is much higher than the minimum value needed to assure compatibility with all post-Newtonian 
experiments, these theories are at present very close to GR. However, their behavior at early times is very 
different from the usual one and can imply a universe expansion rate during the nucleosynthesis process 
several hundred of times faster than in the FRW cosmology. 

The achievement of such pb values can be explained in the following way. The Universe expansion rate 
at the beginning of primordial nucleosynthesis is slower than in GR. This implies a smaller freezing-out 
temperature and, hence, a tendency in these theories to underproduce 4 He. The smaller uj , the larger this 
trend is. The 4 He underproduction can be balanced by considering larger pb values which in principle could 
imply an excessive D burning. But, contrary to all other attempts of modifying the standard hot big bang 
model, the Universe expansion rate obtained in class-4 theories becomes during BBN faster than in GR. 
Consequently, the D burning is not very effective because it occurs in a shortest time and, hence, large 
Pb values are possible. 

Obviously, in order to evaluate quantitatively the light element production in these theories, it was 
necessary to consider a specific form for w($)(Eq. 1). However, this form already reveals a large number 
of models for which Qb is much larger than in GR. Furthermore, it is almost evident that such a result 
should be also obtained for any other scalar-tensor theory, defined by a different coupling function, but 
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implying a similar non-monotonic behavior for the speed-up factor. As a matter of fact, in order to avoid 
both an overproduction of 4 He and an underproduction of D, the only condition is that the speed-up factor 
be smaller than unity at the beginning of BBN and becomes, during the nucleosynthesis process, larger 
than in GR. 

It is important to note, that we do not claim here that the f2s value is necessarily high. The 
essential point that we want to stress is that, in a homogeneous and isotropic universe composed only by 
known particles, primordial nucleosynthesis does not always impose a very stringent bound on the baryon 
contribution to the density parameter. 

Finally, it must be also pointed out that the study of inflation models in the framework of these 
theories, where a form of the coupling function is now known (Eq. 1), and their consequences on the 
primordial density fluctuation spectrum, as well as on the formation of large-scale structures are now 
important open questions which could lead to a new cosmological scenario. 
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Fig. 1. — Space of theories leading to different upper bounds on £1_b: The observational constraints used 
to construct this figure are: 0.221< Y p <0.243, D/H> 1.5 • 1CT 5 , (D+ 3 Hc)< 1.1 • 1CT 4 , 0.7 • 1(T 10 < 
7 Li/H< 3.5 • 10~ 10 , and h G [0.4,1]. The two symbols (square, circle) correspond to the two particular 
theories commented in the text and presented in figures 2. The dashed region represents the space of 
theories with an uppest limit on Q B requiring an oj value smaller than 500. 
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Fig. 2. — Primordial abundances of a) 4 Hc, b) D/H, c) (D+ 3 Hc)/H, and d) 7 Li, as a function of Qb- 
The theories shown in this figure are defined by A 2 = 0.4, logio(A 2 /fc) = 7.8, iv = 1.3 • 10 9 (dashed line, 
corresponding to the circle symbol on figure 1), A 2 = 0.2, logio(A 2 /A:) = 8, ivq = 15.75 • 10 3 (dotted line, 
the square symbol on figure 1) and GR (solid line). A wider observational range for the 7 Li/H abundance 

is also rlisrtl^vprl 
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